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FROM GRAINS TO PLANETESIMALS 

Andrew N. Youdin^ 



Abstract. This pedagogical review covers an unsolved problem in the 
theory of protoplanetary disks: the growth of dust grains into plan- 
etesimals, solids at least a kilometer in size. I summarize timescale 
constraints imposed on planetesimal formation by circumstellar disk 
observations, analysis of meteorites, and aerodynamic radial migra- 
tion. The infall of < meter-sized solids in a hundred years is the 
most stringent constraint. I review proposed mechanisms for plan- 
etesimal formation. CoUisional coagulation models are informed by 
laboratory studies of microgravity collisions. The gravitational collapse 
(or Safronov-Goldreich- Ward) hypothesis involves detailed study of the 
interaction between solid particles and turbulent gas. I cover the ba- 
sics of aerodynamic drag in protoplanetary disks, including radial drift 
and vertical sedimentation. I describe various mechanisms for particle 
concentration in gas disks - including turbulent pressure maxima, drag 
instabilities and long-lived anticylonic vortices. I derive a general result 
for the minimum size for a vortex to trap particles in a sub-Keplerian 
disk. Recent numerical simulations demonstrate that particle clumping 
in turbulent protoplanetary disks can trigger gravitational collapse. I 
discuss several outstanding issues in the field. 



• 1— I 1 Introduction 

X: . 

V^ This chapter serves as a pedagogical introduction to planetesimal fo r matio n and 



1 related processes in protoplanetary disks. I defer to iDominik et al.l ( 20071 ) for a 



comprehensive review of research in this rapidly progressing field. Section [^serves 
as an extended introduction by discussing timescale constraints from disk obser- 
vations, meteoritics, and dynamics. The purposeful overlap with other chapters in 
this volume aims to reinforce the overlap between sub-fields. Section [3] describes 
formation mechanisms: collisional coagulation ( JIS.ip and gravitational collapse 
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( ^3.2p . Section |4] covers the dynamics of gas drag in protoplanetary disks, in- 
cluding the basic form of drag forces f ij4.ip , vertical settling and radial migration 
( ij4.2p . and particle trapping mechanisms (§ 



2 Timescale Constraints 

The agreement between T Tauri disk lifetimes and the inferred formation time of 
meteorites is a remarkable result. Planetesimal formation appears to be a slow or 
ongoing process in the several Myr lifetimes of protoplanetary disks. By contrast 
the radial migration times of cm to m-sized solids can be under a hundred years! 

2.1 Disk Lifetimes 

The standard answer to most questions about disk lifetimes is "a few million years." 
The precise question matters, since main sequence stars (several Gyr old) host low 
luminosity debris disks (see chapters by Beichnian and Kalas) . The Spitzer Spac e 
Telescope detects disks with dust masses as low as ~ 10~^M(p (lMever et al.l . l2008h . 



Erosive collisions between already formed planetesimals supply the observed de- 
bris. 

We are concerned here with the lifetimes of optically thick disks from which 
planetesimals formed. Near-IR excess emission (i.e. above the stellar photosphere) 
from optically revealed young stellar objects (YSOs) indicates a hot inner disk, 
and provides the best statistical determination of protoplanetary disk lifetimes. 
Disk lifetimes are measured by plotting, for many nearby young stellar clusters 
and moving groups, the fraction of stars with IR excess in a cluster against the 
average stellar age for that cluster (estimated by comparing the HR diagram to 
pre- main sequence models). The disk fraction declines uniformly with age from 
near unity to near zero, with a (roughly linear) decay time o f 5 or 10 Myr (for 
the H-K or cooler K-L colors, respectively). See iHillenbrandl ( 2005 ) for details. 



caveats, and references. The fact that we observe a fraction of roughly coeval stars 
with disks indicates that the scatter among disk lifetimes is also several Myr. 

Millimeter and sub- mm emission traces dust in optically thin outer disk regions, 
which provides an estimate of disk mass. The median dust mass in primordial 
disks (around class I and II YSOs which by definition emit excess IR) is (MHnst) ^ 
10~^ Mf7^ « 34M0 ,with some systems ten times as massive (JAndrews and Williams! . 
[20071). For solar composition, the total (gas-dominated) mass extrapolates to 
(Mdisk) ~ 10~^Mq « lOMjup. This is sufficient to form extrasolar giant planets 
(see the chapters by Udry & Eggenberger and Guillot), and is consistent with the 
standard minimum mass solar nebula (MM SN) model, which is itself extrapolate d 
fro m the solid material in our Sola r System (|Weidenschilhng| . ll977tlHavashill98lh . 



Andrews and Williams! ( 20071 ) also found that class III YSOs, which lack IR 



excess, also show little sub-mm emission. Their average dust mass (dominated 
by upper limits) is only (Mdust) < 5 x 10~^Mq k, O.2M0. This supports using 
near-IR excess to measure protoplanetary disk lifetimes. 
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The gas content of disks cannot be weighed as accurately as the dust. Never- 
theless there i s ample evidence o f gas in protoplanetary disks, especially the hot 



inner regions (jNaiita et alj . 120071 ). For instance IR excess correlates strongly with 
accretion diagnostics like Ha emission. We infer that planetesimals form in a gas- 
rich environment so that drag forces must be included in planetesimal formation 
theories. However, there are no strong constraints on dust-to-gas ratios, which 
need not match solar (or extrasolar host star) compositions. 

We conclude that planetesimals should form during the first few Myr of pre- 
main sequence evolution, when disks have sufficient dust to form terrestrial planets 
and cores (see the chapter by Ida). However, the radiative inefficiency of > cm- 
sized solids makes it difficult to observe when planetesimals actually form. 



2.2 Meteoritics 

Solar System (hereafter SS) data provides precise details, unmatched by astronom- 
ical observations, but for a single (admittedly special) system. Primitive meteorites 
- also called chondrites - preserve early SS history. Unlike planets and large aster- 
oids, these rocks were never differentiated into an iron-rich core and rocky mantle. 
They are relatively unaltered since their time of formation. 

Most primitive meteorites contain small (< 1 cm) inclusions which were heated 
to high temperatures (^ 2000 K), thereby removing volatile elements and turning 
amorphous dust crystalline. Spectral signatures of crystalline dust in protoplane- 
tary disks (see chapter by Augereau) show that thermal processing is ubiquitous. 

Chondrules are the most abundant type of inclusion, comprising more than 
70% by volume of the most common primitive meteorites, appropriately called 
"ordinary chondrites." CAIs (calcium aluminum inclusions) are rarer, but are even 
more refractory (i.e. were subject to a more extreme heating event). The forma- 
tion mechanism(s), specifically the heating mechanism, for chondrules and CAIs 
remains a subject of active research. The leading hypotheses are: (1) heating near 
near the young Sun with subsequent outward transport (as in the X-wind model of 
Shu et al.lll996l l. (2) passage though shock waves in the gas disk and (3) co l lisions 



between planetesimals or protoplanets. See the review bv IConnollv et al.l ( 20061 1 



and other chapters in this volume (by Aleon and Wood) for details. 

CAIs are famous for being the oldest objects with precise ages in the SS. The 
record holder is 456 7.11 ± 0.16 Myr old, as determined by 207p|3_206p|3 dating 
( Russell et al.l . l2006l and references therein, hereafter R06). It is said that CAIs 



give the "age of the Solar System," whereas they actually date a specific event (the 
solidification of CAIs) in the history of the Solar System. There are certainly older 
events, s uch as the collapse of the Sun's parent molecular cloud into a protostar 
and disk ( Terebev et al.l . Il984 ). which cannot be dated precisely. There may even 



be older rocks; some iron meteo rites (which come f rom the cores of differentiated 
asteroids) could predate CAIs (JKleine et al.l . |2005[ ). The age difference between 



chondrules and CAIs sets a lower limit to the formation timescale for primitive me- 
teorites, since the inclusions obviously predate their incorporation into the "parent 
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The decay of radioactive ^^Al (ii/2 = 0.73 Myr) provides evidence that CAIs 
are older than most chondrules by a few Myr. Essentially all ^^Al in early solar 
system rocks has now decayed into stable ^^Mg. (Only ~ XO^^^'*" of the original 
^^Al remains after 4.56 Gyr!) Thus short-lived radionuclides (SLRs) like ^^Al are 
termed "extinct" . (The chapter by Williams discusses how SLRs, specifically ^° Fe, 
ii/2 = 1.5 Myr, were introduced to the SS.) The excess of the decay product ^^Mg 
- measured by comparing mineral phases with different Al/Mg ratios - determines 
the amount of ^^Al in the rock "initially," i.e. when it last solidified. Comparing 
the initial (^^Al/^^Al)o ratios between two samples, say a CAI and a chondrule, 
gives the age difference. At, according to exponential decay: 



At = ti/2 ln2 



(26Al/27Al)oa 



(26A1/27A1)„. 



(2.1) 



This technique assumes that ^^Al was uniformly distributed in the disk, otherwise 
differences in (^^Al/^^Al)o could indicate different local formation environment. 
For instance the X-wind model gen erates SLRs near the protostar by spallation 
reactions, not uniformly in the disk ( Gounelle et al.l . 120011 ). 

Most unaltered CAIs share the canonical value (^^Al/^'^Al)o fv i - 5 x 10"^ 
(R06). Applying equation (|2.1|) gives a narrow time span, ~ 0.25 Myr, for most 
CAI formation (if the above caveat on uniform mixing is set aside). By con- 
trast, chondrules with well measured ^^Al (there are only a few) show a range 
of (26Al/27Al)o « 0.24 ± 0.17 X 10"^ — 1.4 ± 0.3 x IQ-'^ (R06). Equation (|2T|) 
implies chondrules formed over an extended period of ~ 2 Myr starting ^ 1 Myr 
after CAI formation. 

Long-lived radionuclides give absolute ages of early SS rocks by using radioac- 
tive isotopes which are not extinct. The most precise ages come from the ^"^Pb- 
^°^Pb technique, which makes use of the decay chains of ^^^U ^^""^Pb (ii/2 = 0.704 
Gyr) and ^s^U ->206pb (i^/j = 4.47 Gyr). Recent measurements give 207pb-206pb 
ages of 4565.7 ± 0.4 Myr, 4564.7 ± 0.6 Myr, and 4562.7 ± 0.5 Myr for chondrules 
in CV, CR and CB meteorite classes, respectively (R06). These best- fit ages are 
respectively 1.4, 2.4 and 4.4 Myr after the oldest CAI mentioned abovecl 

Thus absolute ^"^Pb-^^^Pb dates confirm evidence from relative ^^Al chronom- 
etry that chondrule formation spanned ^3 — 5 Myr after the formation of CAIs. 
R06 summarize the evidence in more detail, including cautionary notes. Thus 
planetesimal formation was either delayed or ongoing for several Myr. This age 
spread fits within the observed lifetimes of protoplanetary disks. Rapid planetesi- 
mal formation in < 1 Myr is not supported by the data. 



^We hope that chondrites are among the first generation of planetesimals formed. This is not 
the case if the coUisional hypothesis for chondrule formation holds. 

^That CAI is from a CV meteorite, so it's not too surprising that it is closest in age to the 
CV chondrule. Note that CV CAIs are large, and thus easier to isotopically date. 
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2.3 Radial Migration by Aerodynamic Drag 

Disk dynamics also constrains the time it takes for planets to form. Planets, and 
their building blocks, have the inconvenient tendency to fall into the star. The 
migration of planets is a significant issue for their survival and orbital evolution, 
especially for hot Jupiters. However the aerodynamic migration of small solids near 
the meter-size barrier imposes the most stringent timescale constraint on planet 
formation: « 100 years. (This is rivaled only by thesis, job and grant deadlines!) 
Moreover the aerodynamic drift mechanism is simple and robust, requiring only 
that the orbital motion of the gas disk is slightly sub-Keplerian due to gas pressure. 
Disk models are generically hotter and denser closer to the YSO, so the requirement 
of a globally decreasing pressure gradient is almost certainly met. 

A solid particle meets a headwind as it orbits in a sub-Keplerian gas disk|j This 
headwind saps the particle of angular momentum causing it to inspiral towards 
the star. Local pressure perturbations (including vortices) can temporarily halt 
the inflow in special locations (see ^4.3p . H owever a net influx will continue with 
a reduced rate (an order unity correction in Ijohansen et al.l . l2006bl ) that depends 
on the efficiency and duration of the particle "traps." 

The amplitude of the sub-Keplerian gas rotation, Sv^fi^g = v^^g — vk (with 
UK = ^jGI^h/ R — QR the Keplerian orbital speed at the cylindrical radius, 
i?), is readily found by balancing the force from the radial pressure gradient, 
fp^R = —p~^dP/dR (where pg and P are the gas density and pressure), with the 
Coriolis force, /cor,_R = 2J7(5w0,g. Setting fp^R + fcor,R = gives: 

1 dP c m 

^'"<P.s = 7, — n^^ = -sp— = -iWK « -56— , (2.2) 

2pgil oR UK s 



where Cg = \J Pj pg is the isothermal sound speed, and r\ is the dimensionless 
measure of pressure support. The factor sp = — (91nP/91ni?)/2 depends only on 
the slope of the gas pressure, with sp w 1.6 in Hayashi's MMSN model. Thus Sv^^g 
is independent of disk mass and depends mainly on the midplane temperature, 
T (X c^. Pressure support is weaker, Svs.^ ~ — 30(i?/AU) ^/^'^ m/s, for passively- 
irradiated flared-disk models (jChiang and Goldreichl . 119971 . hereafter CG97 jf|, but 



stronger in active accretion disks. 

The peak radial drift speed of a particle is simply Sv^^g, directed inwards, 
for particles whose aerodynamic coupling time equals the orbital time (see ij4.2[) . 
Meter-diameter compact solids have the fastest drift speeds at 1 AU in the MMSN, 
while smaller solids drift fastest in the outer disk, e.g. 3 cm at 35 AU. While the 
size of the fastest migrating particle varies with gas density, the peak drift speed 
does not. Aerodynamic drag can bring a particle toward the star in a drift time 

idrift = R/ max(udrift) ~ 85 (i?/AU) yr (2.3) 



•'if the particle orbit is sufficiently eccentric (e > 2r7 — r}'^ to be precise) then it will experience 
a tailwind near perihelion. However eccentricities (and inclinations) not onl y damp more rapidl y 
than semi-major axes, but also increase the orbit-averaged rate of in-spiral IIAdachi et al.L [19761) . 

■*Since CG97 derive T ^ 150(il/AU)-3/7 K vs. T = 280(iJ/AU)-i/2 k for the MMSN. 
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for the MMSN model, or idrift « 160 (i?/AU)"/i4 yr for the CG97 modeli Avoid- 
ing catastrophic loss of planetesimal forming material requires rapid growth though 
the fastest drifting sizes - direct collapse could skip the problematic size range en- 
tirely - and/or efficient trapping of particles (S3)- 

For comparison, I briefly review planetary migration (see also chapters by 
Terquem and Ida). "Type I" migration occurs because a planet launches spiral 
waves in the gas disk. The gravitational back-reaction of the waves on the planet 



exerts tidal torques, which result in a net inw ard migration (Ward, 119971 ). The 



migration timescale from the 3D calculation by iTanaka et al.l ( 



20021) is: 

where Hg — c^j fl is the gas scaleheight. Migration is faster if either the planetary 
mass, Mpi, or the characteristic disk mass, SgR^ (the MMSN is assumed in the 
estimate), increases relative to the stellar mass M^. However, there is a limit. 
Tidal torques can overwhelm disk viscosity and open an annular gap in the disk 



when the planet mass exceeds (iLin and Papaloizoul 1198 



5/8 



Afp,>V40^(f) M* - O.SAW^^ (^J , (2.5) 

where a is the dimensionless parameterization of the disk's turbulent viscosity 
1/ = ac^/f2. Applying inequality (|2.5p to equation (|2.4[) gives an upper limit to 
the type-I migration rate of 



7/8 



Once a planet opens a gap, "type 11" migration occurs on the viscous timescale, 
R^ _ 1 Z' i? ^\,^,, -,„4,.., lO^V ^ \ 

''inis-II '^ ; 



- - -TxlO^yr--^ . (2.7) 



and is generally directed inward (following the flow of the accretion disk). 

This summary glosses over the many uncertainties and difficult technical issues 
associated with planet migration. For instance the horseshoe orbits of gas in the 
corotational region of the planet could lead to a rapid type III migration whose 



direction depends on initial conditions ()Peplihski et al.l 120081 ) . We conclude that 



aerodynamic drift of small solids (sometimes called type-0 migration) is both more 
rapid and based on better understood processes than the migration of planets. 



^The drift or migration timescale, fmig = R/\^r\ (with vji the local radial drift speed) is the 
characteristic time to travel of order the current distance to the star. The time to actually reach 
the star - the integral quantity J„ * dR/v[i{R) — depends on the form of V[i(R). 
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3 Growth Mechanisms 

We now turn to the proposed mechanisms for the formation of planetesimals: (1) 
successive colhsions between sohds which resuh in agglomeration and (2) the grav- 
itational collapse of smaller solids into a larger planetesimal. The two mechanisms 
are not mutually exclusive: dust grains might grow by sticky collisions until solids 
decouple sufficiently from the gas to collapse under their self-gravity. 

3.1 Collisional Coagulation 

Collisions between particles can lead to the rapid growth of plane tesimals (includ 



ing comets in the outer solar system) if particle sticking is efficient (jWeidenschillinei . 



198Ci Il995l ). The key (and difficult) issue is understanding when collisions yield 



particle growth, as a function of particle sizes, impact speeds, orientation and ma- 
terial properties. Significant progress is being made in the laboratory, including 
micro gravity exp eriments in drop towers, sounding rockets and the space shuttle 
(see iBlum and Wurm . 20081 hereafter BW08 for a comprehensive review) . 

Relative velocities between particles regulate collision rates and sticking prob- 
abilities. Turbulent stirring, aerodynamic drift (both radial and azimuthal) and 
vertical settling contribute to collision speeds during planetesimal formation. The 
speeds increase monotonically with particle size, until de coupling occurs at 



meter-sizes, when relative speeds reach ~ 100 m/s. Fig. 4 of iDominik et al.l ()2007l . 
see references therein for details) shows this clearly. 

Collisional growth requires (1) a binding energy (BE) which is significant com- 
pared to the kinetic energy (KE) and (2) a sufficient fraction of KE dissipation, 
/diss, during the impact. Objects should stick when the initial KE obeys 

KE_oo<|BE|-^%^. (3.1) 

1 — /diss 

The logic, illustrated in Fig. [H is as follows. At the moment of impact, the 
KE increases to KE_oo + |BE| by falling into the potential well (it's easy to vi- 
sualize smooth acceleration in a gravitational potential, but contact forces work 
also). Dissipation during the impact then gives KE ^ (1 — /diss)(KE_oo + |BE|). 
Climbing back out of the potential well gives a final kinetic energy KE+oo = 
(1 — /diss)(KE_oo -I- |BE|) — |BE|. Imposing the condition to remain bound, 
KE+oo < 0, reproduces equation (I3.1|l . This heuristic model formalizes the idea 
that sticking requires some combination of gentle collisions, strong binding and 
energy dissipation. Collisions between a small projectile (which dominates the 
center-of-mass kinetic energy) and a large target (which contributes to the po- 
tential energy and provides a sink for dissipated energy) are especially favorable 
for growth. A cruder (but useful) approximation is that coagulation occurs when 
KEi < I BE I, which assumes that /diss is not very close to cither or 1. 

The available binding energies are electrostatic and gravitational interactions. 
Short-range van der Waals interactions explain the sticking of ind ividual grains 
and small aggregates at speeds < 1 m/s (JDominik and Tieleng| . ll997 ). As particles 
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(|) Sticking Bouncing/Fragmentation 

KE 





Fig. 1. A schematic illustration showing the evolution of kinetic (KE) and potential or 
binding energy (BE) during a collision. Sticking requires dissipation of KE so that objects 
remain bound. Escape to infinity in the left example would require negative KE, and 
thus does not occur. Bouncing or fragmentation (right) occurs if there is less dissipation 
and objects escape with positive KE. See equation (|3.f [I and surrounding text. 



increase from /ini to m-sizes, coagulation is more difficu lt since: (1) th e available 
BE declines as the surface area-to- volume ratio drops (JYoudinl . 120041 ) . while si- 
multaneously (2) coUisional KE rises as particles decouple from the gas. Net (and 
opposite) charges can help, but feasibility has not been demonstrated. Experi- 
ments find that mm-sized aggregates do not stick at any speed, with bouncing 
(fragmentation) for impacts above (below) ~ m/s (BW08). 

Gravitational BE is only significant once massive planetesimals have formed. 
The KEi < I BE I criterion is equivalent to requiring (for roughly equal masses) 
that the collision speed is less than the surface escape speed: Wcoii ^ Wcsc ^ 
1 (a/m) mm/s, here normalized to a compact particle with radius a = 1 m. This 
velocity scale is orders of magnitude too small to account for sticking since Ucoii ^ 
10 m/s due to aerodynamic drift alone (assuming a modest dispersion in sizes). 

So what about the scenario of small projectiles sticking to a few "runaway" 
bodies that become planetesimals? The populati on of small bodies coul d be main- 
tained by erosive collisions between large bodies ( Johansen et al.L 120081 ). However 
collisions of small bodies (<C m-sized and tied to the sub-Keplerian gas) with large 
solids (on decoupled Keplerian orbits) face the full brunt of sytematic azimuthal 
drift, Sv^^g > 50 m/s, plus a fluctuating contribution from stirring by turbulent 
gas. Since these speeds are comparable to a sand-blaster (credit to Eugene Chiang 
for this analogy) it is difficult to imagine net growth. BW08 conclude conserva- 
tively that; "The direct formation of kilometer-sized planetesimals cannot (yet?) 
be understood via sticking collisions." This author concurs. 

Many collision mod els predict the growth of bodies > cm in ~ 1000 yr (BW08 
and references therein). iDullemond and Dominitg ( 2005 ) make the important point 
that such rapid growth leads to a premature decline in excess IR and mm emission, 
which would contradict observed Myr disk lifetimes ( ii2.ip . Fragmentation (and/or 
inefficient sticking in the first place) is needed to maintain a population of small 
grains. Astronomical observations are a powerful constraint on our imperfect 
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knowledge of coUisional physics. 

3.2 Gravitational Collapse 

If the "bottom up" growth of planetesimals by coUisional coagulation stalls, then 
a "top down" phase is needed, during which self-gravity fragments the disk of 
small solids into planetesimals. Gravitational instability (GI) is responsible for 
the formation of most objects in the universe, including cosmological structure 
and stars. It has been hypothesisized that giant plane ts can form by gravitational 
fragmentation of protoplanetary gas disks ( Bossl . ll99'iT). However the gas is unlikely 



to cool sufficiently rapidly to collapse into planets (JRafikovl . 120071 ). Here we are 
concerned with the collapse of the solids in protoplanetary disks, which cool (i.e. 
damp their random kinetic energy) rapidly by gas drag and inelastic collisions. 

3.2.1 Stability Criteria 

Safronovl (|l969f ) and iGoldreich and WardI (|l973l . hereafter GW73) independently 



proposed that planetesimals form by GI. They applied two stan dard crit e ria for 
gravitational collapse: the Roche (or critical) density and the iToomrd ( 19641) 



Q-parameter. The Roche limit is usually invoked for tidal disruption, e.g. the 
rings inside Saturn's Roche lobe, but collapse is merely the inverse of disrup- 
tion. Gravitational collapse is possible when the self-gravity, with an acceleration 
/sg ~ Gp£ for a region of density p and size £, exceeds the differential tidal gravity, 
/tide ~ GM^,£/R^, at a distance R from the central mass, M*. This gives the simple 
collapse criterion p > pr ~ M^/R^, but no information about the characteristic 
lengthscale since we ignored pressure stabilization of small scales. 

The Toomre criterion states that thin disks are gravitationally unstable when 

Qt ^ ^ < 1 (3.2) 

where c is the random velocity (for a ga s the soun d speed), fi is the orbital fre- 
quencjo, and E is the surface density (see lShul . ll992l or lBinnev and Tremainell2008 



for a derivation). This criterion says that mass promotes gravitational collapse, 
while pressure and epicyclic restoring forces (represented by c and fi respectively) 
stabilize the disk. The dispersion relation (see above texts) restricts the unstable 
wavelengths to be less than At = Att'^GS / ff' , since angular momentum conserva- 
tion prevents large regions from collapsing. The Roche criteria can be rearranged 
to look like the Toomre criteria: 

where the scaleheight, H k, S / p. Equations p.2p and (|3.3p are roughly equivalent 
for c sa ilH , which holds for gas disks and most p article orbits, but does not h old 



for small particles suspended in the gas disk (e.g. lYoudin and Lithwickl . 120071 ) . 



'More generally the epicyclic frequency, k, should appear, but k = i? in a Keplerian disk. 
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The standard GI criteria (eqs. |3.2| and |3.3| ) apply to single componen t disks . 



In th e presence of gas drag, collapse is possible even when Qt and Qr 3> 1 (jWard . 



2000l though the basic idea is in GW73)13 The transfer of (excess) angular mo- 



mentum, from the particle layer to the gas, allows longer wavelengths to collapse. 
Since self-gravity is weaker, these dissipative gravitational instabilities can be slow, 
especially for small particles that must push through the gas. 

3.2.2 Turbulence: A Curse and a Cure 

The stirring of small solids by turbulent gas is the fundamental obstacle to grav- 
itational collapse. In the absence of any turbulence, particles would settle to 
arbitrarily high densities in the midplane - faster than they drift radially - so that 
GI is trivial. Hoping that turbulence disappears, even just in the midplane, is not 
a very satisfactory solution. For instance, turbulent stresses can propagate int o 



"dead zones" from the active surface layers of accretion disks (|Oishi et al.l . 120071 ). 
Moreover, particle settling toward the midplane robustly generates midplane 
turbulence. As solids sediment into a dense layer, they drag the midplane gas along 
at faster (closer to Keplerian) speeds, generating vertical shear with the overlying 
sub-Keplerian gas. This shear in turn triggers Kelvin-Helmholz instabilities (KHI) 
which develop into turbulence. The presence of KHI was recognized by GW73 
and noted to be an obstacle to GI by Weidenschillind ( 19801 ). Research on KHI in 



protoplanetary disks continues (e.g. IChianei |2008[) 



Particles can also generate midplane turbulence via streaming instabilities. 



^articles can also generate miapiane turouicnce via streaming mstaDiiities, 
hereafter SI (JYoudin and Goodmanl . l2005[ ). Streaming motions between solids and 



gas - the naturally occuring radial and azimuthal drift ■- supply the free energy 
for SI. By contrast gas velocity gradients, i.e. shear, are responsible for KHI. To 
model SI, the gas and solid components must evolve indep endently, either as two 
fluids , or with a hybrid hydrodynamics and N-b odv scheme IjYoudin a n d JohanserJ . 



2007^. A hvbrid of K HI and SI (simulated by Ijohansen et al.l . l2006al in 2D and 



^ 1 1 3 

Johansen et al.l . l2007l hereafter J07, in 3D) operates in real disks with stratification 



and imperfect particle-gas coupling (i.e. with slippage allowed). 

Midplane turbulence may not be the fatal flaw for GI, for two reasons. Both 
effects are augmented in disks with high metallicity (and thus large initial solid 
abundances), consistent with the observation that exoplanets are more common 
around metal-rich stars (Fischer and Valentil. 2005). Fi rst, KHI can only stir flnite 



amounts of sohds (jSekival . ll998tlYoudin and Shul . l2002l ). If the local abundance of 



solids relative to gas exceeds ~ 10 times solar abundances, then any excess solids 
cannot be stirred by vertical shear, and will be subject to GI (absent sufficient 
stirring by other sources of turbulence). This saturation limit exists because the 
amount of destabilizing vertical shear is limited (by Sv^,g « 50 m/s) while the 
stabilizing anti-buoyancy increases as more solids are poured into the midplane. 



^These models do include radial turbulent diffusion, which slows but does not prevent this 
mode of collapse ( Youdin, in prep) . 
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Second, while turbulence diffuses particles in an average sense, it can also gen- 
erate intermittent particle clumps. Enhanced self-gravity in an overdense clump 
is conducive to GI. Mechanisms for particle concentration in turbulence (broadly 
defined to include longer- Hved vortices) are described in ^4.31 Particle clumping is 



particularly strong for SI (jJohansen and Youdinl . 120071 ). since particle overdensities 



generate perturbations to the drag force that sustain the growth of SI. 

The models of J07 demonstrate the viability of planetesimal formation by GI 
in turbulent disks. Their hybrid MHD/particle-mesh simulations include several 
sources of turbulence - vertical stratification (for KHI), particle-gas streaming 
(for SI), and magnetic fields to drive the magneto-rotational instability (MRI)|f| 
Particle clumps grew in the absence of self-gravity, and then collapsed rapidly when 
self-gravity was included. The largest bound clump accreted at an astonishing rate 
of ~ 0.5Mceres pGr Orbital period. A clump mass that exceeds a nominal km-sized 
planetesimal by orders of magnitude is not a serious concern, since hierarchical 
fragmentation in the later stages of collapse (not yet modeled) is plausible. 

The initial particle sizes in J07 are in the 15 — 60 cm range, making them 
moderately decoupled from the gas. Growth by coagulation to these sizes is not 
certain ( §3.ip . Strengthening the GI (or Safronov-Goldreich-Ward) hypothesis 
requires (1) a better understanding of particle growth beyond ~ 10 cm and/or (2) 
pushing the dynamical models to smaller initial particle sizes, which should give 
slower growth rates and may require higher initial particle abundances. 

4 Aerodynamic Coupling in Particle-Gas Disks 

4.1 Drag Laws and the Stopping Time 

While planets interact with gas disks gravitationally, smaller solid bodies are cou- 
pled to the gas disk primarily by drag forces. The drag force, Fd, on a particle of 
mass rUp moving with a speed Av = Vp — Vg relative to the gas is 

Fd = -mpAv/tstop, (4.1) 

where istop is the damping timescale for particle motion relative to the gas. The 
form of istop depends on particle properties - such as the internal density, ps, and 
spherical radius, a - and on properties of the gas disk - pg and Cg - as: 

iX = /'sa/(PgCs) ifa<9A/4 (4.2a) 

4*°pk- ^ i^top ■ 4a/(9A) if9A/4<a<A/(4Ma) (4.2b) 

^Pp • {ajXfl^ Ma-^/V4 if A/(4Ma) < a < 200A/Ma (4.2c) 

t*-^ = t'^l^ ■ 6/Ma if a > 200A/Ma (4.2d) 



^stop 



*The supplement to J07 shows that GI occurs in the absence of magnetic fields, though higher 
particle surface densities were needed. This is surprising, since (naively) less turbulence should 
favor GI. The clumping in the MRI turbulence strengthens the feedback of drag forces on the 
gas, which in turn enhances the concentration by SI. 
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Fig. 2. Dimensionless stopping time, Ts, vs. disk radius in the midplane of a minimum 
mass nebula model for particle sizes from fira - km. The sharp transitions between 
different drag regimes are approximations. 

where Ma = |At;|/cs, A ex l/pg is the gas mean free path, and Re = 4aMa/A wiU 
be the Reynolds number of the flow around the particle. The cases are written in 
order of increasing particle size: Epstein's Law of drag from molecular collisions, 
Stokes' Law for viscous drag when Re < 1, an approximate intermediate Re case, 
and the drag from a fully developed turbulent wake for Re > 800. 

We measure the dynamical significance of drag forces with the parameter. 



= nt 



stop ■ 



(4.3) 



For Ts ^ 1 particles are carried along with the gas, while for Tg 3> 1 gas drag is a 
small correction to Keplerian orbits. Fig. [2] plots Tg values for several particle sizes 
in the midplane of the MMSN with a gas surface density, Sg = 2000 {R/ AV)-^/'^ 
g/cm , and midplane temperature, T = 280(i?/AU)^^/^ K. The relative veloci- 
ties needed for the turbulent drag formulae (eqs. |4.2cl4.2d] ) are determined self- 
consistently by inplane drift (eq. [4.7| ) . Turbulent stirring and vertical oscillations 
will give a small correction. 



4.2 Settling and Drift 

4.2.1 In-plane Drift 

We now derive the particle drift motions. First the orbital speed of the gas, v^f^^g, 
is given by centrifugal balance with the Keplerian gravity and pressure forces. 



JR = -vl/R - p-^dP/dR . 



(4.4) 



We ignore the back reaction of drag forces on the gas, valid in the test-particle limit 
when the density of particles pp ^ pg. When pressure forces are small compared 
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Fig. 3. (Left:) Particle radial drift timescale for several sizes vs. disk radius. The fastest 
migration at any radius is for a particle with Ts = 1 {grey curve). (Right:) Particle settling 
times, which are longer than the orbital time (grey curve), but shorter than the radial 
drift times. 

to radial gravitjothe approximate solution, Sv^^g « i'0,g — fK = —V^k ~ —56 m/s, 
matches equation (|2.2p . Rotation is again sub-Keplerian for the (usual) outwardly 
decreasing pressure, with 77 = —{dP/dR)/(2pgf2'^R) « (cs/wk)^ ^ 10~^. 

The particle response satisfies the radial and azimuthal equations of motion: 



R-R(t)^ = 
R'(j) + 27?0 = 



—v^/R — R/tstop 
-{R<j)- W0,g)/istop 



(4.5a) 
(4.5b) 



where the radial gas velocity vanishes. For steady state solutions, we set R ~ 
and let eee f2{R) + 6v^^p/R so that R(f> « -(3/2)127?. Using equation ([^ for 
the gas velocity gives the particle drift speeds. 



■^fl.,] 



R = - 



1+t2 



(4.6) 



and (Su^.j 



'??i'k/(1 



The inward drift speed is maximized at — 77UK for 



Ts — I particles. The total (in-plane) particle speed relative to gas. 



Av\ ^ d R-^ + {Sv^^, 



VVk) 



l + r2 



rjVK 



(4.7) 



asymptotes to \Av\ = tjvk for large particles (ts ^ 1) on Keplerian orbits. 

Fig. [3](/e/t panel) plots the radial drift timescale, tdrift = -R/|i'fl,p| for a range 
of particle sizes. The fastest drift, tdrift, min = iv^)^^ = 88(i?/AU) yr is for Tg — 1 
particles {shown in grey). 



^Otherwise the disk would be spherical! 
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4.2.2 Vertical Settling 

Vertical particle motion decouples from in-plane motion, obeying the simple ODE, 

Z = -(i-W2,g)/istop-i^^2:, (4.8) 

where g^ ~ —fi'^z is the Keplerian vertical gravity near the midplane [z <C R). 
When forcing by the vertical gas velocity, w^.g, vanishes (for a laminar disk) we 
have a damped harmonic oscillator, with a damping - i.e. particle settling - time: 

Re[l - y/\ - 4r2] ^t^ 

The final, approximate expression is a simple matching of the tight and loose 
coupling limits. For r^ <C 1, overdamped oscillations bring particles directly to 
the midplane at their terminal velocity Wsctt(-z) = gztstop = ~^'^ztstop, giving a 
settling time tgctt = z/vsctt{z) — (l^^tgtop)^^- For large particles with Tg :» 1, 
oscillation amplitudes are damped as exp[— i/(2istop)]- Fig.^ {right panel) plots 
settling times for several particle sizes. Note that fim-sized grains have igctt ~ Myr, 
of order the gas disk lifetime, emphasizing the need for some particle coagulation. 
Forcing by the turbulent gas maintains a finite particle layer thickness, Hp. 
This is given by balancing tgctt with the turbulent diffusion timescale across Hp: 
idiff ^ H^/D y We skip details of how the particle diffusion coefficient. Dp, 



depends on Ts (jYoudin and Li thwictc. ' 20071) . The particle scaleheight is given by 



a single simple expression (jCarballido et al.l . 120061 : lYoudin and Lithwickll20 07i). 



Hp 



(^^f^s, l) H, , (4.10) 



where a^ = Dg/{H^n) parametrizes turbulent diffusion in the gasl^ The ^/a^J^^ 
factor tells us that weaker turbulence and/or larger, heavier particles result in a 
thinner particle later. Taking the minimum simply prevents Hp > Hg. 

4.2.3 Small Particles Move (Relatively) to High Pressure 

A simple rule of thumb is that particles move to regions of high pressure. Midplane 
settling and inward radial dr ift (or falling rair jj here on Earth) follow this pattern. 



To derive this result (see also lJohansen et al. 



2006bl ). and get a simple prescription 



for the drift speeds of small particles, consider the equations of motion for the gas 
velocity, Vg, and a particle labeled by index i: 

^-fU^) = -— (4.11a) 

Dt pg 

t-stop t-stop 



^"^This Og is not identical - but can be similar in magnitude - to the usual a for angular 
momentum transport, which includes Maxwell stresses. 
^^ Sadly there was little snow during the Winter School. 
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Whenever Tg <C 1, particles are frozen- in to the gas to lowest order, so that 
Xi « Vg{xi) + 0{ts). Thus the left hand sides of eq. [4.11b ..b] are in balance, 
including external forces, fext - such as gravity or Coriolis forces - which affect 
both components. Thus the "extra" force felt by the gas - pressure gradients - 
must balance the "extra" force felt by the particles - drago Consequently 

VP 

^v,^ i3top + 0(T,2), (4.12) 

Ps 

which gives the correct radial drift and vertical settling speeds, but not the order 
Tg azimuthal drift speed. We apply this result to vortex trapping (eq. [4.16J ). 

4.3 Particle Trapping Mechanisms 

A promising route to planetesimal formation involves the clumping of particles by 
gas structures in the disk, which are mostly turbulent in origin. Particle density 
enhancements can then seed gravitational collapse (i j3.2p . These trapping mecha- 
nisms also slow the migration of solids, a serious problem in its own right (i i2.3|) . 

4.3.1 Pressure Reversal 

The tendency for particles to drift radially inward is again due to the azimuthal 
headwind caused by the global radial pressure gradient, dPo/dR ^ —pgC^/R < 
(the subscript "o" now indicates a smoothing of fluctuations on scales f <C i?). 
A stationary point, with no headwind or radial inflow, exists wherever a local 
pressure perturbation, P' , flattens the total radial pressure gradient, 

dP/dR\,^,^, = dPo/dR + dP'/dx = , (4.13) 

with X the local radial coordinate. Pressure maxima (not minima) are stable 
stationary points. Small pressure perturbations, P' ^ Po, reverse the gradient if 
they vary on length scales, i < RP' / Pq <C i?. A lower limit, (. > 4top = Rotstop = 
2r]RT^/{l + Tg), ensures that particles can stop at the stationary point. Long-lived 
fluctuations are more effective traps, with igtop the minimum duration. The fastest 
drifters (with Tg = 1) concentrate most efficiently when these criteria are met. 

What structures can generate the requ ired pressure ma xima? Self-gravitating 
gaseous spiral arms are one possibility (e.g. iRice et al.l . l2004l) . The pressure pertur- 
bations must be substantial for large (AU-scale) density waves (applying the gra- 
dient argument above). This scenario requires a massive, and presumably young 
disk, which may not be favored according to the arguments in (f|2|). 

Pressure fluctuations in disk turbulence are another candidate. Simple scalings 
give the pressure perturbations in terms of the turbulent Mach number and a as: 



^^Pressure (also, magnetic fields, ignored here) has a negligible effect on dust grains and solids 
because of their large "molecular weight." Drag forces have little effect on the gas in the test 
particle limit when the mass density of particles Pp <S^ pg- 
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P' /Po ~ Mat^ri-, = w^^j.,^/cg ~ a, for eddies with a dominant length scale (. ~ ^/oiHg 
and turnover time ioddy ~ f2~^. Both the pressure reversal and particle stopping 
(applied to Ts < 1 since larger solids will fly through these eddies) criteria give the 
same requirement for trapping: a > (Hg/R)'^ ~ 10"'^, a fairly standard value. 

The effectiveness depends on details beyond these simple estimates. Inverse 
cascades to long-lived, shear-elongated struct ures will enhanc e trapping. Some 
MRI simulations have found strong clumping (jjohansen et al.L l2006bl . 120071 ). and 



work is ongoing to understand the robustness of these results. 

4.3.2 Vortex Trapping 

The possibility of particle trapping in long-lived vortices receives considerable at- 
tention, at least in part due to the beautiful analog y with Jupiter's red spot. W e 
ignore the relevant question of vortex stability (see iBarranco and Marcud . 120051 ) . 



and derive two results: (1) the trapping of particles in anti-cylonic vortices and 
(2) a constraint on vortex size in the (realistic) case that the gas rotation is sub- 
Keplerian. The latter result may be new. I will not disc uss concentration b y vor- 



tices at the dissipation scale of Komogorov turbulence (ICuzzi et al.l.l200lh. sin ce 
the time and spatial scales may be too short to be relevant ( Cuzzi et al.l . l2008h . 

Consider a region of constant z-vorticity; lu — dxUy — dyllx, where x and y are 
the local radial and azimuthal coordinates in a frame rotating uniformly with the 
vortex center. The velocity field in the vortex must satisfy 

Ux = enyy , Uy^-^x, (4.14) 

where e is the x:y axis ratio, and i?v is the frequency of fluid oscillations. The 
vorticity is w = —i7y{l+€'^)/e, but the the excess vorticity we = oj—luk = (jJ+3f2/2 
is the more relevant quantity in (here Keplerian) shear flows. 

We me ntion two specif ic vortex solutions, but will keep our derivations gen- 
eral. The ISaffmanI ( 19921 and references therein) vortex solution with ily = 



(3/2)/2e/(l — e) matches onto the Kepler shear flow Uy = —(3/2)fix for i nfinite 



azimuthal elongation: e — > 0. The "planet" solution of iGoodman et al.l (|l987r i 
with i7v = ef2 1/3/(1 — e^) exactly solves the compressible continuity equation. 

To find the particle trapping criterion we start with the pressure forces in the 
vortex. Balancing against tidal gravity, Coriolis and incrtial accelerations gives 



1 dP n dUr- / , 2i7i7v 9 

_^— = 3n^x + 2nUy-Uy—^^ (sn^ - + n^ 

Pg ox oy \ e 



1 f)p fjTJ 

— — = ^2nUx-Ux^ = {~2enQy + n^)y. (4.15b) 

Pg oy ox 

The divergence of the particle velocity, Vp, then follows from equation (|4.12p 

V • Vp = V • At> = ^fi^V^P = 2ts {lue + ^v/^) • (4-16) 
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Thus particle concentration requires we < —fl^jfl < 0, i.e. anticylonic relative 
vorticity. Our result only holds for Tg <C 1, but even for the loose coupling only 
anticylonic vortices trap particles. See lChavanid ( 2OO0 Hbr details. 

the centers are not sta- 
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For vortices embedded in sub-Keplerian gas disksL 
tionary points for particle trapping. No drag forces would be exerted there to 
enforce a sub-Keplerian particle orbit. The answer - that the stationary point 
remains radially centered on the vortex, but moves forward in azimuth - is ob- 
vious in retrospect. In the Keplerian frame, particles trapped in a sub-Keplerian 
vortex experience a radial Coriolis force, fcoi,x = —2f2rjVK- This must be bal- 
anced by the radial drag force of the vortex flowing over the stationary particle as 
fcoT,x + Uooiys)/tstop = 0. The stationary point is 

and the final approximate expression is for an elongated (e ^ 1) Saffman vortex 
(i.e. a weak perturbation to the Keplerian vorticity). The constraint that j/s lie 
inside the vortex imposes a minimum size for particle trapping, which becomes 
more stringent for longer stopping times. 

4.3.3 Drag Instabilities 

A surprising consequence of drag forces in disks is that they drive instabilities 
which clump particles. This possibility was first explored bv lGoodman and Pindor 



(|200d hereafter GPOO). Instead of calculating microscopic drag forces as in equa- 



tion ()4.ip . they developed a single fluid model for a (height-integrated) particle 
layer subject to turbulent boundary layer drag. They found a robust clumping 
instability for various formulations of turbulent momentum transport. The cru- 
cial factor - which they demonstrate with an intuitive "toy model" - is that the 
frictional drag acceleration must vary with the surface density of the particle layer. 



tionai drag acceleration must vary witn tne suriace density ot tne particle layer. 
The streaming instability (SI) of lYoudin and GoodmanI pOOSl l is also quite 



robust and leads to clumping, but the setup of the model is different. The calcula- 
tion is a bit more detailed, including the 3D coupled motions of both particles (a 
pressureless fluid in the analytic treatment) and gas. However the instability arises 
spontaneously from a laminar background, and the required physical ingredients 
are minimal: (1) local Keplerian dynamics (axisymmetric and unstratified in the 
simplest form), (2) drag forces which act on both particles and (by Newton's Third 
Law) gas, and (3) a global pressure gradient. The microscopic drag acceleration on 
a particle, f^^ = F^/mp = — Au/tstop is independent of the spatial mass density 
of particles, pp, unlike the turbulent boundary layer drag of GPOO. Instead it is the 
back reaction on the gas, /^ g = — (Pp/Pg)/d,p = (Pp/Pg)^''^/^stop, which allows 
particle clumping to mediate the instability. Linear perturbations to pp interact 
with the background aerodynamic drift, Av, which was described in ii4.2.1l 



^''Local vortex solutions are unaffected by a constant global pressure gradient, except for 
Galilean velocity transform by — jjrix. 
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The ultimate energy source for all drag instabilities - the secular instability of 
GPOO, SI and KHI - is pressure-supported, sub-Keplerian gas motion. While the 
above discussion emphasizes robustness, drag forces do not replace the need for 
an angular momentum transport mechanism in accretion disks (MRI turbulence is 
currently the leading candidate). Drag affects the gas only where mass loading is 
significant, pp > pg, i.e. in the sedimented midplane. Also drag effects are strongest 
for moderate coupling (near Tg ~ 1). Instabilities are weak for small particles 
"slaved" to the gas, or for large bodies that mostly ignore the gas. Finally, at least 
for SI, an gular momentum transport i s inwards, the wrong direct ion for stellar 
accretion (JYoudin and Johansenl . l2007t I Johansen and Youdinl . 120071 ). However, as 



discussed in ^3.2.21 drag instabilities play a crucial role in planetesimal formation 
theories, especially if clumping leads to gravitational collapse. 

5 Concluding Remarks 

Planetesimal formation remains unsolved, despite considerable progress on several 
fronts. Laboratory studies of collisions place constraints on sticking assumptions 
( Blum and Wurml . [2008). Numerical simulations coupled with analytic theory 



provide insights to physical processes in disks. The response of particles to - and 
the effect of particles on - turbulent gas is an especially promising area. Contrary 
to previous expectations, simulations show that gravitational coll apse of > 10 cm 



solids into planetesimals occurs in turbulent protoplanetary disks (jJohansen et al 



20071 ). Forming planetesimals directly by gravitational collapse from chondrule- 
sized, < 1 mm, solids would be an elegant solution. But such small particles are 
more tightly coupled to the gas, making it difhcult to understand (analytically or 
numerically) how they might collapse in standard disk models. 

Sadly numerical (and even laboratory) experiments lack direct observational 
confirmation. But the initial conditions and basic ingredients of theoretical models 
can be compared to the meteoritic record and observations of protoplanetary disks. 
Spitzer is delivering exciting, detailed results. Soon Herschel and ALMA will 
further expand our knowledge of the conditions in - and structure of - disks. 
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